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Summary. - We study a R 2 model of gravity with torsion in a closed Friedmann- 
Robertson- Walker universe. The model is cast in Hamiltonian form subtracting from 
the original Lagrangian the total time derivative of JkIr, where fx is proportional to the 
trace of the extrinsic curvature tensor, and fa is obtained differentiating the Lagrangian 
with respect to the highest derivative. Torsion is found to lead to a primary constraint 
linear in the momenta and a secondary constraint quadratic in the momenta, and the full 
field equations are finally worked out in detail. Problems to be studied for further research 
are the solution of these equations and the quantization of the model. One could then try 
to study a new class of quantum cosmological models with torsion. 
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1. - Introduction. 

A central issue in modern cosmology is the role played by theories of gravitation other 
than Einstein's general relativity in getting a better understanding of the early universe. 
Among the various alternative theories, the ones with torsion [1] are still receiving careful 
consideration. Most theories with torsion are such that spin can be thought of as the 
source of torsion. There are at least two very important motivations for studying this kind 
of gravitational theories with torsion (we are indebted to C. Stornaiolo for explaining these 
issues). In fact: 

a) The holonomy theorems [2,3] imply that torsion and curvature are related, respec- 
tively, to the groups of translations and of homogeneous transformations in the tangent 
vector spaces to a manifold. In general relativity the spin does not play any dynamical 
role, so that the infinitesimal holonomy group is given just by the homogeneous trans- 
formations. The introduction of torsion related to spin introduces a much stronger link 
between gravitation and particle physics, because it extends the holonomy group to the 
translations. A very enlightening discussion of gauge translations can be found in [4] (to- 
gether with the literature given therein) and in [5]. In particular, the introduction of [5] 
makes clear from the very beginning the main geometric role played by the translations 
in the gauge group: they change a principal fibre bundle having no special relationship 
between the points on the fibres and the base manifold into the bundle of linear frames of 
the base manifold. A different view is expressed in sect. 5 of [6]. 

b) At the very high densities present in the early universe, the effects of spin are no 
longer negligible. In addition, theoretical investigations have shown that torsion may lead 
to the avoidance of the big-bang singularity in some cases (see for example [7]). 

A good treatment of a theory where the sources of gravity are mass and spin is for 
example [8]. In that paper the Poincare gauge theory of gravitation is studied, setting up a 
first-order Lagrangian formalism in a Riemann-Cartan space-time. However, in our paper 
we will not focus our attention on this kind of theories. In fact, we are more interested in 
getting a preliminary understanding of some basic features of the canonical structure of 
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R + R 2 theories of gravity with torsion. The reader should keep in mind that in theories 
with torsion derivatives higher than 2 are not present if frame and connection are regarded 
as independent variables. The R 2 terms arise already in the torsion- free case as quantum 
corrections to the effective action. In general such theories suffer from serious problems due 
to the presence of ghosts. However in [9-11] it has been found that some R + R 2 theories 
are ghost-free. This is indeed a delicate point. In fact in studying teleparallel theories 
(where the constraint of vanishing curvature is imposed), the authors in [12] found that 
the most general quadratic Lagrangian for such theories leads to dipole ghosts. These are 
field modes with free-field energies not bounded from below, and exist in view of p 4 -poles. 
The authors in [12] also studied the general Poincare gauge theory, and in sect. 5 of their 
paper they discussed in detail the difference between their work and the approach used in 
[9] in dealing with p 4 -poles. Other reasons for studying R + R 2 theories with torsion are: 

i) The calculation of the Hamiltonian for R + R 2 models with torsion is a very inter- 
esting field of research, because it provides an extremely fertile interplay between Dirac's 
theory of constrained Hamiltonian systems and the problems of gravitation and cosmol- 
ogy. Here we would like to mention the important paper [13], where Dirac's method is 
applied to the most general, parity-conserving Lagrangian, which is at most quadratic in 
the torsion and the curvature (the special R 2 case is treated in [14], whereas higher-order 
theories in the Riemannian case are studied for example in [15]). 

ii) In the torsion- free case, R + R 2 theories have been shown to lead to an inflationary 
universe [16,17]. Moreover, other theories with torsion can lead to inflation (see for example 
[18,19]). 

iii) In R + R 2 theories in the most general case, torsion may be expected to propagate, 
a remarkable property which, however, is not shared by the Einstein-Cartan-Sciama-Kibble 
theory based on spin. 

In sect. 2 we study a R 2 model of gravity with torsion in a closed Friedmann- 
Robertson- Walker (hereafter referred to as FRW) universe. The consideration of this 
model is suggested by a Lagrangian studied in [10]. We cast the model in Hamiltonian 
form using a canonical variable obtained differentiating the Lagrangian with respect to the 
highest derivative. Torsion makes itself manifest in that it induces one primary constraint 
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01, linear in the momenta. Requiring the preservation in time of 0i, we end up with a 
secondary constraint 4>2 quadratic in the momenta. Contact is also made with the canoni- 
cal structure of the corresponding torsion-free model. In sect. 3 we summarize the results 
obtained, and we mention the problems to be studied for further research. 

2. - R 2 theories with torsion: a closed FRW model. 

We are here interested in a R 2 theory of gravity with torsion. The Lagrangian of our 
model is assumed to be given by 167tGL g = N\fh fi^^Rj , where no assumption on the 
sign of n is made. Indeed, a slightly more general Lagrangian would be 



This model was at first studied by Rauch and Nieh [10], though not in a FRW universe and 
not using canonical techniques. Thus our Lagrangian is a special case of the Lagrangian 
studied in [10] when A = 0. In so doing our analytical derivations will be simplified, which 
in turn implies an advantage in trying to solve numerically the field equations, which 
will be shown to be rather involved. Moreover, setting A = we realize a more direct 
comparison with the canonical structure of the torsion-free model studied in sect. 3 of 
[17]. Thus both our formalism and the physical content of our model are rather different 
from what has been discussed in [10]. The only basic ideas we strictly need to recall are 
the following. In a theory of gravity with torsion we have a four-dimensional space-time 
manifold with a metric tensor and a nonsymmetric linear connection T x ^ v which obeys the 
metricity condition [1]. The torsion tensor in a coordinate frame is defined by 





(2.1) 



One also defines a contorsion tensor 



< A c A c A 



+ 



(2.2) 
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so that 



[IV 



-C 



(2.3) 



where j £ v j are the Christoffel symbols. We shall use a convention according to which the 
curvature tensor is given by [20] 



pa pa pa I pA pa _ pA pa 



(2.4) 



We will focus our attention on a closed FRW universe. In view of the hypothesis of 
spatial homogeneity and isotropy, we are interested in a torsion tensor whose only 
nonvanishing components are 



(2.5) 



In our model, the metric may be locally cast in the form 



ds 2 = -N 2 (t)dt 2 + a 2 (t) (d X 2 + (sin X ) W + (sin X ) 2 (sin#) 2 # 2 ) . (2.6) 



Thus one finds (see appendix) 



= + 12 ( — 
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(2.7) 



We now define 



16ttG 



/ (^R) 2 Ny/hd l x 



Therefore, putting f3 = 37r/2G, a = lg(a), n = /j/3 : we find 



/ 



L dt 



(2.8) 



(2.9) 



L = ^ exp[3a]iV- 



6exp[— 2a] + 



12a 2 



aQ 96Q 2 6 d (a-2Q 



(2.10) 
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Thus, putting 

Q y, rfNdt, (2.11) 



N 

denoting by a prime the derivative with respect to r, and defining fx = 3/i, we obtain the 
relations 

/ = Jldr, (2.12) 

~ r 9 i 2 

L = fi° exp[3a] exp[-2a] + 2a' - \2ya! + 16y 2 + a" - 2y' . (2.13) 

Now, in the torsion-free case (compare with [17]) one defines the variable z = dL/da", 
and one considers the Lagrangian (we are indebted to D. Giulini for clarifying this point): 

L' = L- 4-(<*'z) ■ (2-14) 

The geometrical meaning of z is that it is proportional to the scalar curvature through 
fx exp[3a]. In defining (2.14), we fix a and z rather than a and a' at the initial and final 
times. In our model with torsion, we can no longer use (2.14) as our final Lagrangian V . 
In fact, if one defines p a = dL'/da', p y = dL'/dy' and p z = dL'/dz' using (2.14), the 
resulting relations do not involve y', whereas p y = —2z. Thus it is not possible to compute 
the Legendre transform: H = p a a' + p z z' + p y y' — L'. However, we can point out that in 
(2.14) a' is a function proportional to the trace of the extrinsic curvature tensor. In our 
model with torsion, a' gets replaced by a' — 2y (see (2.7) and (2.13)). This suggests to 
define 

d 



(a'-2y)z . (2.15) 



In deriving (2.15) in a more systematic way, we can point out that the total derivative 
appearing in (2.14) is 

„ dL , d dL 
da" dr da" 
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Thus, if torsion has to play a role, we expect having to add the term 



,dL 



d dL 




dr dy' 



so that we end up with definition (2.15). To sum up, in defining the correct L' we require 
that 

a) it must be possible to compute the Legendre transform: H = p a a' +p z z' +p y y' — L' '; 

b) V must reduce to (2.14) in the torsion-free case, when y = 0; 

c) L' must not contain derivatives of y higher than y', because these derivatives are 
absent in (2.13); 

d) L' must be defined in a unique way, through a general method; 

e) V must be of the kind V = L — (d/dr)(fKfR) where /# is the function a' — 2y, 
proportional to the trace K of the full extrinsic curvature tensor, and fn is proportional 
through Jl exp[3a] to the full scalar curvature ^R. 

It is also worth remarking that: 

i) in view of (2.1)-(2.3) and (2.5), the only components of the connection (2.3) which 
depend both on a and on Q are Tq 1: Tq 2 and Tq 3 , and they are all equal to a — 2Q = 



iii) if we require that the constraint p y ~ should be avoided and that all constraints 
produced by torsion should be linear in the momenta, we are led to define 



N(a'-2y); 

ii) formula (2.15) can be written in the form 




which is formally of the same kind of (2.14), if we define x = a 



2fydr; 




But in so doing there are infinitely many forms of L' , one for each choice of the set of 
functions {fi(a)}, and our construction is completely arbitrary, violating requirement d). 
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Thus we will use (2.15), where, putting 

V = u' , 



(2.16) 



one has 



dL n 



exp[-2a] + 2a' 2 - Ylv! a! + 16u' 2 + a" - 2u" 



(2.17) 



As one can easily check, putting y = u' one has the advantage of dealing with less involved 
calculations. Using the definitions p a = dL' /da 1 ', p z = dL'/dz', p u = dL'/du', we find 
that the effective Hamiltonian H is given by 

# = # + 701. (2.18) 

In (2.18), one has 

H = p a a' + p z z' + p u u - L' 

, 2 . VzPu . z 2 exp[-3a] 
= -4zp z + ~2~ + ^exp[-2a] (2.19) 

with the primary constraint 

(/) 1 =p u + 2p a - 4zp z . (2.20) 

The constraint (pi must now be preserved in time using Dirac's method [21] for 
constrained Hamiltonian systems. Namely, computing the Poisson bracket of 0i with H, 
we find the secondary constraint: 

7z 2 

4> 2 = 16zp 2 z - 2p u p z + — — g exp[-3a] - 8zexp[-2a\ . (2.21) 

LI 

In other words, torsion introduces a primary constraint 0i linear in the momenta, which 
in turn leads to one more secondary constraint quadratic in the momenta. Finally, so 
as to preserve in time the constraint 4> 2 ~ 0, we require that the following Poisson bracket 
must vanish: 

{^ 2 ,#} = {(/> 2 ,tf} + 7{02,<M • (2.22) 
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Setting the right-hand side of (2.22) equal to zero and solving for 7 we find 

{02, 



7 = 



{02, 0l} 



where 



2 49^ 2 



{02, 0i} = 64zp z - 8p u p z q~~ exp[-3a] + 64zexp[-2a] , 



{fo,H}= ^96exp[-2 a ]^- 72eX P[ M z^jp. 
+ (^ 9eX ^ 3a ^ - 6exp[-2a]^ p u . 



z 



Using the constraint 0i 0, we can also cast (2.19) in the form 



H = -2zp 2 z - p a p z + - e ^| Q - ^exp[-2a] 



Thus the equations of motion of our model are given by 



/ 9H 

a =dp- a = ~ Pz + 2 ^ 



1 dH 

z = —— = -4zp z - p a - A^z , 

Op z 

, dH 

u = j— = 7 , 
op u 

dH 3z 2 r ■. . n . 

Pa = ~~da = 4^0 exp ^ ^ ~~ 2^exp[-2aJ , 

. dH n 9 zexpl — 3a] r „ . 

P* = — ^ = IpI + exp[-2a] + A 1Pz , 



dH 
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plus the constraints: (pi « 0, (p2 ~ and the formula for 7 (see (2.19)-(2.21) and (2.23)- 
(2.25)). In other words, we deal with a coupled system of six first-order ordinary differential 
equations subject to constraints. A remarkable result is that we have written the field 
equations in Hamiltonian form without having to solve any differential equation involving 
torsion. Moreover, (2.32) implies that p u = const. It is also very important remarking 
that, imposing the particular value H — (see (2.19)), we can solve for zp 2 . Inserting this 
relation into the constraint <p2 = 0, we obtain the relation 

exp[a] - ^z = . (2.33) 

In other words, requiring the constraints (pi = 0, (p2 = and setting H = is equivalent to 
require that $1, $2 and H must vanish, where $1 = (pi and $2 is given by the left-hand 
side of (2.33). The constraints (pi and (p 2 are second class (they cannot be brought into 
the first class by suitable linear combinations [21]). It is also worth emphasizing that in 
the torsion-free case no primary constraints such as (pi are found to arise. Setting y = 
in (2.13) and using (2.14) one finds then 

a' = -Vz , (2.34) 

z' = -4zp z - p a , (2.35) 
3z 2 

Pa = ex P[- 3 «] - 2^exp[-2a] , (2.36) 

p' z = 2p 2 z -^ exp[-3a] + exp[-2a] , (2.37) 
plus the particular value for the Hamiltonian here chosen 

z 2 

2zp 2 z + p a p z + z exp[-2a] - — ^ exp[-3a] = . (2.38) 

So as to avoid confusion, it is important to remark that, using our notation, one has 

, a „ 1 d ( a \ 

a = No: ' a = NJt liv^y ' 
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whereas in other papers such as [17] one defines 

, da da dt 
a = — — = — — — = aa = a , 
dr) dt drj 

.. da' da 1 
a = — — = a—— = aa . 
dr) dt 

System (2.27)-(2.32) with the constraints (2.20)-(2.21) is very difficult, and at the classical 
level we have not a criterion which picks out a preferred choice of initial conditions. In- 
deed, exact solutions for theories with torsion are known for a large class of theories. An 
important recent paper is [22], where the authors have studied the Hamiltonian structure 
of Poincare gauge theory. In the case of spherical symmetry and for the charged Taub- 
NUT metric, the authors of [22] have obtained the most general torsion configuration for a 
large class of quadratic Lagrangians (see also references in [22] for related work). However, 
our original aim was just the study of the canonical structure of the model described by 
(2.5)-(2.6) and (2.8) at the classical level. Thus we prefer to stop here our analysis, and 
to outline in detail the unsolved issues in the next section. 

3. - Conclusions. 

In this paper we have studied a Hamiltonian approach to a closed FRW universe with 
torsion. The consideration of our R 2 Lagrangian (see (2.8)) has been suggested by a model 
at first studied by Rauch and Nieh [10]. Our main results are the following: 

1) the technique used by Horowitz [17] in the torsion- free case can be generalized using 
(2.15)-(2.17); 

2) torsion modifies the canonical structure of the theory in that it induces a pri- 
mary constraint linear in the momenta and a secondary constraint 4>2 quadratic in the 
momenta (see (2.20)-(2.21)); 

3) the field equations are much more complicated, and are given by (2.27)-(2.32). 
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Thus the most remarkable feature of our model seems to be the appearance of the 
additional constraints <f>i and 2 , or equivalently $1 = </>i and $ 2 (see (2.33)). In gen- 
eral relativity, the secondary constraints reflect the invariance of the theory under four- 
dimensional diffeomorphisms. Thus in our model we can say that torsion introduces an 
additional invariance under four-diffeomorphisms, leading to 02, which contributes to the 
generation of the dynamics through 7 (see (2.23)-(2.25) and (2.27)-(2.32)). In the as yet 
undeveloped minisuperspace model, we expect the physical state ^ will still be a function 
of the three-metric hij and of the extrinsic curvature Kij as in the torsion-free case [17]. 
In fact, the nonvanishing part of torsion is in our case just the one contained in (see 
appendix of [20]), and from (2.27) and (2.29) one has a' — 2u' = a' — 2y = —p z , which is 
formally identical to (2.34), because in our model with torsion the trace K is proportional 
to a' — 2y (see sect. 2). Indeed, Hamiltonian methods had also been used for the canonical 
quantization of more complicated R 2 theories of gravity in the torsion- free case [23,24], 
and in studying spherically symmetric gravitational fields with electromagnetism and a 
massless scalar field as sources [25]. For theories with torsion, a rather important paper 
is [26], and in fact the ansatz (2.5) of our paper is closely related to the ansatz studied 
in sect. 2 and sect. 6 of [26]. A minor difference between our work and [26] lies in the 
different parametrization of the FRW metric, and in the fact the authors of [26] set the 
lapse function equal to one, whereas we keep it arbitrary. But the main difference is the 
following: the authors of [26] do not perform a canonical analysis, and thus they do not 
deal with six first-order equations subject to constraints. A comparison with their analyti- 
cal derivation has not yet been possible because it is very difficult to solve our equations of 
motion. Therefore, the main problems to be studied for further research are the choice of 
the initial conditions for the numerical integration of the field equations, the quantization 
of the model using path integrals and the attempt of applying the Hartle-Hawking proposal 
[27] or other proposals. In fact we cannot solve the field equations without knowing the 
initial conditions. Quantum cosmology can provide a way for doing so, but then one has to 
work out the path-integral quantization, discussing in detail the measure in the path inte- 
gral, the metrics we are summing over, the gauge fixing and the ghost action (namely the 
action involving Faddeev-Popov ghosts). Moreover, it is necessary to understand whether 
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the underlying quantum theory is affected by negative-norm states, whose occurrence can 
be a serious drawback. The quantization procedure is here complicated by the additional 
constraints $1 and $2- 

Indeed, torsion is much less ad hoc than the fundamental scalar fields considered so 
far in cosmological models (for the role played by scalar fields, see for example [28,29]), 
though its existence is still to be proved. Moreover, in the torsion- free case R + R 2 theories 
may play an important role not only for the inflationary era, or in giving a self-consistent 
probabilistic interpretation of the wave function of the universe [30] , but also in completely 
different processes such as black-hole evaporation [31]. This is why we think it is important 
to try to solve the above-mentioned problems. To be honest, we also must say that it 
is not yet clear to which extent we can use a closed FRW metric so as to study the 
early universe. Therefore the model we have studied is just a mathematical idealization. 
Nevertheless, the derivations obtained show that it is worth studying theories of gravity 
with torsion in further detail from a Hamiltonian point of view. We hope that our paper can 
stimulate further interaction between canonical gravity, theories with torsion and quantum 
cosmology. We think it would be also very useful to study the singularity problem for R 2 
theories with torsion both in the isotropic and in the anisotropic case. 

* * * 

We are grateful to P. D'Eath, D. Giulini and J. Louko for useful conversations and to 
St. John's College for financial support. 

APPENDIX 

We are interested in the formula for the scalar curvature ^R for a theory of gravity 
with torsion. This kind of calculation is performed in detail in [32]. However, this reference 



13 



HA MIL TONIA N STRUCTURE OF A FRW UNIVERSE WITH TORSION 



is not easily available. Thus we prefer to rely on the appendix of [20]. In that appendix, 
Pilati works out the Lagrangian density ignoring a total divergence. Now, in general form, 
and bearing in mind eqs. (A. 16)- (A. 18) of [20], we can write 



(4) 



R = K {ab) K^ - K [ab] KM + K 2 — 2r (ab) K^ + 2r [ab] K^ + 2Kr - 2 
+ i3) R + A + B + C, 



K 

N 



(AA) 



where 



A = h 13 



V^rai ),j ~ K^ji ),m "I" <^ji <^ml ~ ^mi^jl 



+ 



I 



(A.2) 



B = ~^N\ ij -N i K^ , 



C = -p a q a + 2p a la -2p a p a 



(A3) 
(AA) 



In (A.4), p a and q a are the ones defined in (A.30)-(A.31) of [20]. Using (2.5)-(2.6) and 
(A.1)-(A.4) of our appendix, one finds (2.7). Other useful references for the Hamiltonian 
formulation of ECSK theory are [33-35]. Finally, we wish to remark that, putting 



S = 



2h kl S km 
N 



(A.l) leads to a surface term involving torsion in the action integral 



8tvG 



x 



dM 



in addition to the usual one of Einstein's theory [36-38]. 
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